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Abstract 

Three-dimensional neutral hydrogen mapping using the redshifted 21 cm line has recently 
emerged as a promising cosmological probe. Within the framework of slow-roll reconstruction, 
we analyze how well the inflationary potential can be reconstructed by combining data from 
21 cm experiments and cosmic microwave background data from the Planck satellite. We con- 
sider inflationary models classified according to the amplitude of their tensor component, and 
show that 21 cm measurements can significantly improve constraints on the slow-roll parameters 
and determine the shape of the inflationary potential. 
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1 Introduction 



Inflation is the prominent paradigm of the early universe that explains the flatness over cosmological 
scales, the Gaussianity of density perturbations and the near scale invariance of the cosmic power 
spectrum. Accelerated cosmic expansion during inflation pushes perturbation modes from casually 
connected scales to outside the horizon. After re-entering the horizon these superhorizon modes 
provide homogeneity over apparently casually disconnected scales, and give rise to the peaks in the 
power spectrum of the cosmic microwave background (CMB) radiation which has been measured 
with unprecedented precision over a five-year period by the Wilkinson Microwave Anisotropy Probe 
(WMAP5) pQ. The Planck [2 J satellite, planned to be launched in 2009, and continued observation 
by WMAP will further exploit the rich information from both CMB temperature and polarization 
power spectra. 

However, the mechanism that drives the early universe into inflation remains an open question. 
Generically inflation can be modelled by an inflationary field rolling down a potential 0E]. Models 
may be large field [5], small field [6] and hybrid [Tj and have been widely studied. Alternatively 
an inverse method [8] focuses solely on the kinematics of rolling and reconstructs the inflationary 
potential in a model-independent manner. The slow-roll parameters are defined in terms of the 
derivatives of the potential. These parameters can determine the primordial power spectrum that 
sheds light on how well the inflationary potential can be experimentally probed. Slow-roll param- 
eters have been utilized lately to analyze inflation with WMAP data [9] and the upcoming Planck 
project [TOlfTT]. 

A number of radio telescopes are currently being proposed, planned or constructed to observe the 
redshifted 21 cm hydrogen line from the Epoch of Reionization (EoR), e.g., MWA [T2], 21CMA [T5] , 
LOFAR [14J, GMRT [IS], PAPER pi], Square Kilometer Array (SKA) [17], and Fast Fourier 
Transform Telescope (FFTT) [18] . 21 cm tomography maps the neutral hydrogen in the universe 
over a wide range of redshifts and provides a promising cosmological probe, with arguably greater 
potential than CMB and galaxy surveys. Several studies have investigated the accuracies with which 
cosmological parameters can be measured by upcoming 21 cm experiments, both by mapping diffuse 
hydrogen before and during the EoR [19] and by mapping neutral hydrogen in the galactic halo 
after reionization [20]. In particular, the FFTT experiment optimized for 21 cm tomography can 
improve measurement of the cosmological parameters to an unprecedented level [21] . Consequently, 
precision measurements from 21 cm tomography open a new window to constrain inflation in the 
early universe. 

In this paper, we adopt a model-independent approach and forecast how accurately the shape 
of the inflationary potential can be reconstructed by combining the 21 cm data from FFTT or SKA 
and the CMB data from Planck. In the next two sections we outline the reconstruction method 
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and assumptions about the 21 cm power spectrum. In Sections H] and we describe the two classes 
of kinematical models and their analysis. We display our results in Section [H 



2 Potential reconstruction 

We briefly outline the potential reconstruction method and refer the reader to Refs. [231 [23] for 
extensive discussions. 

Consider a flat universe whose energy-momentum tensor is dominated by an inflaton field (j) 
evolving monotonically with time in a potential V{<f>). With the Hubble parameter H expressed in 
terms of <f>, the equation of motion of 4> and the Friedmann equation can be written as 



mpi 

4-7T 



H'W, (1) 



and 



i - 



(2) 



where m-p\ is the Planck mass, primes and overdots denote derivatives with respect to (ft and time 
respectively, and 

ml x (H'{cj>)\ 2 



4tt \H(<P) 
Inflation occurs so long as e < 1. 

A series of higher order parameters are obtained by successive differentiation [25] : 

where n > 1 and the usual slow-roll parameters are rj = X\ and £ = A2. No assumption of slow-roll 
is made in the definition of these parameters. If the hierarchy of differential equations is truncated 
so that X n = for n > m, an exact solution for H((f>) (up to a normalization factor) can be 
found [26j. Once H((f>) is known, the shape of the potential V(</>) is determined. 

The evolution of the slow-roll parameters is conveniently expressed as a function of the number 
of e-folds before the end of inflation N. With 

dN Y- ^ (5) 
the flow of the slow-roll parameters is given by [3| 

-j^ = [(n - l)Ai - ne] X n + A n+ i . (7) 
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To solve these equations, we need to specify values of the slow-roll parameters when observable 
modes left the horizon. We denote these by a "0" subscript and take ko = 0.05 Mpc -1 to be the 
fiducial mode. We set (j>Q = 0. 

The spectral indices and their running that define the commonly used power-law parameteri- 
zation of the primordial scalar and tensor power spectra [27\ 

k \ n s -l+|«ln ^ 



P,(k) = A^-j , (8, 

f h \ nt 

p,(k) = A t y , (9 ) 

can be related to the slow-roll parameters. To second order, expressions for parameters that will 
be relevant to our study, are [23] 

n s = l + 2r/o-4eo-2(l + C)eg-i(3-5C)eo% + ^(3-C)Co, (10) 



dn s 1 dn s 

a 



dink l-e dN 

where C = 4(ln2 + 7) — 5, with 7 ~ 0.577, and the tensor to scalar ratio r = A t /A s is 



(11) 



r = 16e [l + 2(-2 + In 2 + 7 )(e - %)] . (12) 

WMAP5 data support a red-tilted (n s < 1) spectrum and r < 0.25 pp. With A s fixed by observa- 
tion, the parameter r determines the tensor amplitude. If r > 0.1, tensor modes are detectable by 
Planck [2H]. 

3 21 cm power spectrum 

We briefly describe the essential background of 21cm cosmology in this section, and refer the 
interested reader to a comprehensive review in Ref. [29J. The redshifted 21 cm line due to the 
neutral hydrogen hyperfine transition can be measured in terms of the brightness temperature 
relative to the CMB temperature [30] . 

rp 1 x _ 3c 3 hA 10 n H (x)[T s (x) -T C mb] n „x 
bW 32<irk B v$T s ( X )(l + z) 2 dv ll /dr ' 1 j 

where Aiq is the spontaneous decay rate of 21 cm transition, nn is the number density of the 
neutral hydrogen gas, T$ is the spin temperature and dv\\/dr is the physical velocity gradient along 
the line of sight (with r the comoving distance). The temperature fluctuation can be parametrized 
in terms of the fluctuation in the ionized fraction 6 X , the matter density fluctuation 8, and the 
gradient of peculiar velocity along the line of sight dv r /dr. During the EoR, the hydrogen gas is 
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heated well above the CMB temperature [31] . so that in the approximation T s 3> TcmBj 

n = p^[i- (xiKi + s x m +s)(i- ^ d ^) , (14) 

\xh) V -" a " r / 

where Xi = 1 — xh is the ionized fraction of hydrogen gas and xh is the fraction of neutral hydrogen. 
The total 21 cm power spectrum P A r(k) is defined by (AT*(k)AT b (k')} = (2vr) 3 (5 3 (k - k')PAT(k), 
where AXfc(k) is the deviation from the mean brightness temperature and k is the comoving wave- 
vector that is the Fourier dual of the real coordinate position r. We restrict our considerations to 
linear perturbation theory (<5 <C 1) and write the Fourier transformed spectrum to leading order as 

^AT(k) = P (k) + P 2 (k)i? + P 4 (k)fi 4 , (15) 

where the multipole coefficients can be written as 



Here fj, = kn is the cosine of angle between the wave- vector and the line of sight. The power spectra 

of matter and ionization fluctuations are denoted by Vss = f^(xH) 2 Ps5, PxS = T^(xH){xi)Ps x 5, 

and V xx = f£(xi) 2 Ps x s x , where f b = (tJ^Tcmb ^h) w l^/T' We account for ionization effects by 
parameterizing the ionization power spectra as [21] 

V xx (k) = b 2 xx [1 + a xx {k R xx ) + (k R xx f] _2fl Vss , (19) 
V x8 (k) = b 2 xS exp[-a x s(kR x s)-(kR x s) 2 ]Vss, (20) 

where b xx and b 2 xS are the amplitudes of the spectra, R xx and R x $ are the effective sizes of the 
ionized bubbles (HII regions), and a xx , ^ xx and a x s are spectral indices. We adopt the fiducial 
values of Table III in Ref. (21]. 



4 Model classification 



P = V ss -2V xS + V; 
P2 = 2{Vs5-V xS ), 
Pa = Vss- 



(16) 
(17) 
(18) 



Kinematically different potentials can be categorized based on the relative sizes of slow-roll param- 
eters. The parameter e plays a critical role that determines the duration of inflation, the rate of 
change of <fi, how much the inflationary potential V((f)) rolls down from its initial height, and the 
tensor to scalar ratio. We follow a recent classification that is based on the size of e pTJ. 
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4.1 High e models. 

High e models yield r > 0.1 so that tensor modes are detectable by Planck. 

One-parameter models, e is the sole parameter in these models and determines the primor- 
dial spectra. As the only free parameter, e is stringently constrained by 21 cm and CMB data. 
However, this single parameter scenario is not easily realized in particle physics. 

Two-parameter models. In these models rj contributes to the evolution equations. Two- 
parameter models resemble a ACDM cosmology with significant tensor power. 

Three-parameter models. Inflationary rolling is described by e, rj and £. These models 
resemble a ACDM model with measurable r and a large £ can contribute significantly to the running 
of scalar spectral index a, breaking scale invariance of the power spectrum. The non-zero £ allows 
the rolling to speed up at late times and gives a variety of shapes for the potential. £ contributes 
significantly to a when it is numerically comparable to the other two parameters. Generically, £ 
speeds up the evolution of e and a large £ causes a prompt end to inflation with small N. 

4.2 Low e models. 

In these models e is vanishingly small when ko leaves the horizon. We set eo = 1CT 8 . This represents 
extremely slow rolling at horizon-crossing. In such models non-zero higher order parameters cause e 
to grow super-exponentially near the end of inflation and the potential falls abruptly with a cliff-like 
feature. 

Two-parameter models. These models resemble ACDM with near scale invariance in the 
power spectrum and negligible tensors. The parameter rj can be strongly constrained but the 
number of e- folds are generally large because an efficient accelerating mechanism is absent. Within 
95% C. L. constraints from WMAP5, we find that these models give N > 180. A large N indicates 
that inflation must end via a hybrid transition. 

Three-parameter models. A non-zero £ parameter speeds up rolling, significantly lowers 
the number of e-folds and allows a non-zero a. These models can easily be distinguished from 
the two-parameter case. It is noteworthy that in these models it is possible for rolling to be even 
slower than in two-parameter models during most of the inflationary period. This is followed by 
significant late-time acceleration which causes the overall effect of £ to be a speed-up of rolling. 
The phase of slow evolution also occurs in models with higher order kinematical parameters. 

Here we do not investigate low e models with higher order parameters (A n with n > 3) since 
such models are indistinguishable from the three-parameter model. 
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5 Analysis 

21 cm experiments do not directly measure k or Pat s ). The power spectrum Pat( u ) is evaluated 
in the observer's pixel u that is the Fourier dual of the observed vector = 6 x x + 9 y y + Avz where 
(Ox, @y) gives the angular location on the sky plane, Av is the frequency difference from the central 
redshift of a data bin and the z-axis is along the line of sight. By using Pat( u ) instead of Pat s ), 
we avoid the Alcock-Paczynski effect [32J, which arises from the model dependence in the projection 
of the physical wave-vector k over cosmological distances. 

We employ the Fisher matrix formalism to determine the precision of parameter estimation. 
Following Ref. [21], we resolve the 21 cm spectrum Pat( u ) into pixels and the 21 cm Fisher matrix 
is constructed as 



where 5Pat( u ) is the power spectrum measurement error in a pixel at u and A is the combined 
set of cosmological and ionization parameters. 

We consider 21 cm measurements in the redshift range 6.8 — 8.2 with three redshift bins centered 
at z = 7.0, 7.5 and 8.0, with a nonlinear cut-off scale k max = 2 Mpc -1 , and 16000 observation hours. 
Non-Gaussianity of ionization signals is assumed to be negligible in our analysis. We assume that 
the foreground can be perfectly cleaned above the scale k m [ n = 2n/yB where yB is the comoving 
line-of-sight distance width of a single redshift bin. This assumption was shown to be a good 
approximation in Ref. |19j . We consider two detector arrays, SKA and FFTT, which have optimal 
signal-to-noise ratios among planned 21 cm experiments. We assume an azimuthally symmetric 
distribution of baselines in both arrays. The design of SKA has not been finalized. We adopt the 
"smaller antennae" version of SKA, in which the array will have 7000 10 m antennae. We assume 
that 16% of the antennae are concentrated in a nucleus within which the area coverage fraction is 
close to 100%; 4% of the antennae have a coverage density that falls as the inverse square of the 
radius; and 30% are in the annulus where the coverage density is low but rather uniform out to a 
5 km radius. We ignore the measurements from the sparse distribution of the remaining 50% of the 
antenna panels that are outside the annulus. FFTT is a future square kilometer array optimized 
for 21 cm tomography as described in Ref. [T5]. Unlike other interferometers, which add in phase 
the dipoles in each panel or station, FFTT can obtain more information by correlating all of its 
dipoles. We assume that FFTT contains a million 1 m x 1 m dipole antennae in a contiguous core 
subtending a square kilometer, and providing a field-of-view of 2ir steradians. 

The Fisher matrix formalism for the CMB is well established [33] ; for Planck data we follow the 
latest experimental specifications [2j. We include both temperature and polarization measurements 
and assume l ma x = 3000 with three frequency channels while the other channels are used for 




(21) 
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foreground subtraction. The CMB power spectra's parameter dependence is computed using the 
Code for Anisotropics in the Microwave Background (CAMB) |34j . 

The Fisher matrix is cosmology dependent and we work in the fiat (fl^ = 0) standard ACDM 
model and fix f^/i 2 = 0.0074 (neutrino density) and Y p = 0.24 (helium abundance). The fiducial 
values of the non-slow-roll parameters are set near the best-fit of the WMAP5 result [lj: h = 0.72 
(Hubble parameter Hq = WOh km s _1 Mpc -1 ), r = 0.087 (reionization optical depth), Q\ = 0.742 
(dark energy density), f^/i 2 = 0.02273 (physical baryon density), and A s = 0.9. We fix Vss(k) in 
Eqs. ()19p and (|20p when varying cosmological parameters, so that constraints arise only from the 
Vss terms in Pq, P2 and P4. 

The Fisher matrices depend on A that includes (n s , r, a) in Y planck and (n s , a) in F 21cm 
in addition to the non-inflationary parameters. We marginalize over the latter to obtain ^f^l-a) 
and F 2 ^ cr ^. The Jacobian matrix d\ spe c/ d\ sr (where the subscript "spec" indicates (n s , r, a) 
for Planck and (n s , a) for 21 cm experiments), can be used to obtain the Fisher matrices for the 
slow-roll parameter set \ sr = (e, ry,^), 



d\ sr J tspec d\ s ■ 



The three independent spectral parameters allow the Jacobian matrix a maximal rank of three, 
and directly constrain up to three slow-roll parameters. We consider Planck and 21 cm data 
independently, so the combined Fisher matrix is the sum of the contributions, 



-plot pizlcm _|_ -pflancK fOQ^ 

r sr ~ r sr ~r r sr ) V Zt) 7 

which we use to construct a x 2 function, 

X 2 (X sr ) = SlF^S Xsr , (24) 
where d denotes the deviations from the fiducial values of the slow-roll parameters. 



6 Results 

We forecast constraints on the slow-roll parameters at the fiducial points of Table [1] that are 
consistent with WMAP5 results. To supplement the uncertainties listed in the table, we provide 
the corresponding (approximate) uncertainties for the more familiar spectral parameters. The 
joint SKA+Planck (FFTT+Planck) analysis gives the la uncertainties 5n s = 0.0031, 5a = 0.0032 
(5n s = 6 x 10~ 4 , 5a = 2.7 x 10" 4 ). These results roughly apply to all the classes of models in 
Section^ These uncertainties are larger, but consistent with those in Ref. [21] since we marginalize 
over all other parameters, while in Ref. [21 j . r and a are held fixed in computing uncertainties for 
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Table 1: Uncertainties on slow-roll parameters for models classified according to the size of e. The fiducial values 
at the time of horizon-crossing are chosen to be consistent with the 2a ranges favored by WMAP5 data [T]. 
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Figure 1: 2a forecasts for the fiducial points in Table [T] from a Fisher matrix analysis of SKA+Planck (dashed) and 
FFTT+Planck (solid) in the (n s ,r) and (n B ,a) planes. 



n s , and r is fixed in computing uncertainties for a. For r large enough to be measured by Planck, 
^ r |r~o.i = 0.022 and if r is tiny, a bound Sr\ rr ^Q = 0.005 can be placed at la. 21 cm data do not add 
any information on tensor modes. Figure [1] shows this information pictorially. The fiducial points 
are chosen so that the allowed regions have minimal overlap. A comparison of the constraints from 
the joint analyses with that from Planck data alone is shown in Fig. [2j The constraining power 
of 21 cm data comes from their sensitivity to n s and particularly a. 21 cm and Planck data are 
complementarity in their sensitivity to a and r. 

It should be mentioned that higher order corrections to the brightness temperature power 
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0.940 0.945 0.950 0.955 0.960 0.965 
« s 

Figure 2: The impact of 21 cm experiments on parameter estimation. 2<r regions from an analysis of Planck alone, 
21 cm alone, and 21 cm+Planck. The fiducial point for the high e two-parameter model is used. FFTT and Planck 
are complementary: FFTT has good sensitivity to a but no sensitivity to r and Planck has good sensitivity to r but 
not a. We do not show the FFTT+Planck ellipse since it is indistinguishable from the ellipse for FFTT alone. 

spectrum (Eqs. fT5l — I20j) may lead to errors as large as 0(1) in the power spectrum at small scales 
k > 1/iMpc -1 when the neutral fraction is (xu) ~ 0.5 [22]. However, since interferometer array 
measurements are more sensitive to small k modes than to large k modes because of thermal noise, 
cosmological constraints depend only weakly on k max , the nonlinear cutoff scale above which we 
ignore 21 cm contributions to cosmology. Figure 6 of Ref. [21] shows that in their setup, the 
uncertainty in the tilt measured by the FFTT and the Planck data varies from roughly 0.0003 to 
0.0006 to 0.0009 as k max is reduced from 2 Mpc" 1 to 1 Mpc" 1 to 0.6 Mpc -1 . Regardless of the 
exact value of A; max that can be determined by further careful modeling, it is qualitatively robust 
that cosmological constraints from FFTT and Planck data will reach unprecedented precision, e.g., 
the measurement of n s at the level of 10~ 4 . 

To implement Monte Carlo reconstruction of the potential, we randomize slow-roll parameters 
inside the 2a regions allowed by 21 cm+Planck data as the values when the scale k$ left the horizon. 
We then evolve Eqs. ([6]) and (|7|) forward in time. Those cases are selected in which inflation ends 
with the number of e- folds N that pass a prior N m i n < N < N max . The prior on N is necessary 
because (i) a sufficiently large ./V is required to be consistent with the observed horizon size; (ii) 
a small N indicates relatively fast rolling which suggests that higher-order parameters may not 
be small enough to be truncated; (iii) a large N indicates that rolling is extremely slow so that a 
hybrid mechanism might be responsible for end the inflation. While our framework supposes that 
observable inflation is dominated by a single scalar field, it does not preclude the possibility of a 
hybrid transition caused by other fields ending inflation. We use two priors, 40 < N < 70 and 
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Figure 3: Bands of reconstructed potentials at 2a from SKA+Planck (left) and FFTT+Planck (right) for two sets of 
priors on the number of e-folds, 40 < N < 70 (upper) and 30 < N < 500 (lower). Note that the low e two-parameter 
model requires N > 180 and is eliminated by the 40 < N < 70 prior. The unshaded bands in the upper panels are 
shown only for comparison. It is remarkable that the reconstruction using FFTT+Planck is barely affected by the 
c-fold prior. The enlarged SKA+Planck band for the low e three-parameter model for 30 < iV < 500 is a consequence 
of £ ~ being allowed at 2a. 



30 < N < 500. The first prior is typical for a plausible expansion history of our universe with 
Ref. [35J arguing for N between 50 and 60. This first prior does not account for a hybrid transition. 
Our second prior is rather conservative 30 < N < 500, with the large values suggesting that some 
other mechanism brings an abrupt end to inflation. 

In Fig. [3l the bands show the envelopes of possible potentials at the 2a C. L. for each class of 
models with fiducial values as in Table[TJ The envelopes capture the shapes of the potentials because 
the reconstructed potentials do not show any fine dependence on (j). While the reconstruction from 
SKA+Planck is clearly affected by the e-fold prior, it is striking that the reconstruction from 
FFTT+Planck is essentially unaffected. The low e two-parameter model is inconsistent with the 



11 




Figure 4: The impact of 21 cm experiments on potential reconstruction. 2<r bands from Planck alone and 
21 cm+Planck. The left panel is a magnified view of the low e three-parameter model; the SKA+Planck band 
is almost identical to that for Planck alone and is not shown. 



40 < N < 70 prior since WMAP5 data yield N > 180 for these models. The SKA+Planck band 
for the low e three- parameter model expands greatly for the 30 < N < 500 prior because £ ~ 
becomes allowed at 2a. Note that detection of tensors by Planck is not sufficient to guarantee 
satisfactory potential reconstruction using Planck data alone. For example, FFTT data crucially 
improve the reconstruction of the high e 1 parameter model. In Fig. HI we compare results for 
models which require 3 slow-roll parameters for their description. FFTT data narrow down the 2a 
bands considerably. 

We conclude by emphasizing that a joint analysis of 21 cm measurements from FFTT with 
Planck data will significantly help pin down the slow-roll parameters and determine the shape of 
the inflationary potential. The improvement over the reconstruction using Planck data alone may 
stimulate major developments in our understanding of the particle physics responsible for inflation. 
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